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ABSTRACT 

We investigate the redshift evolution of the molecular gas mass fraction (/moi = 
A/j+A/ 1 ’ ■''^here Mmol is molecular gas mass and M* is stellar mass) of galaxies in 
the redshift range of 0 < 2 < 2 as a function of the stellar mass by combining CO 
literature data. We observe a stellar-mass dependence of the /moi evolution where 
massive galaxies have largely depleted their irrolecular gas at 2 = 1, whereas the /moi 
value of less massive galaxies drastically decreases from 2 = 1. We compare the ob¬ 
served M* — /mol relation with theoretical predictioirs from cosmological hydrodynamic 
simulations and semi-analytical models for galaxy formation. Although the theoreti¬ 
cal studies approximately reproduce the observed mass dependence of /moi evolution, 
they tend to underestimate the /moi values, particularly of less massive (< 10^° M 0 ) 
and massive galaxies (> 10^^ Mq) when compared with the observational values. Our 
result suggests the importance of the feedback models which suppress the star for¬ 
mation while simultaneously preserving the molecular gas in order to reproduce the 
observed M* — /moi relation. 
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1 INTRODUCTION 

Current galaxy-formation models are constructed so as to 
reproduce the observed stellar mass function (SMF) of 
galaxies. The statistical properties of the stellar compo¬ 
nents of galaxies have been investigated by means of large 
surveys in the optical and near-infrare d wavelengths (e.g. 
ICole b ^ I 2 OO 1 I : iKochanek et al.ll2001^ . It is known that 
the galaxy formation models based on the canonical cold 
dark matter hypothesis (the so-called ACDM hypothesis) 
overpredict the number densities of less massive and mas¬ 
sive galaxies when compared wi th the observed SMFs of 
galaxies (e.g. iBenson et al.11^031'1 . To resolve this discrep¬ 
ancy, several feedback processes have been proposed to blow¬ 
out and/or heat up the cold-gas co mponent, which is the 
raw material for star fo rmation le.g. Isilk fc Mamonll20l3 : 
ISomerville fc Davd 1201^ 1. For low-mass galaxies, the feed¬ 
back processes of t he blowing-out of the cold gas by su¬ 
pernovae (SNe; e.g. lLars^ll974l : iDekel fc Silklll986ll . pho- 
toionisation by cosmic ultraviolet background radiation (e.g . 
Ilkeuchil 1 19861 : lEfstathi^ 1 19921 : iNavarro fc Steinmet j 119971 : 
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lGnedirJl2000l : [Okamoto et al.ll2008l l. and preheating of the 
inter galactic medium (IGM) by g ravitational pancaking 
le.g. lMo et alll2005l: ILu et al.l[20lii have been considered. 


while for high -mass galaxies, acti ve galactic nucleus (AGN) 
feedback (e.g. Croton_eLal 200^) and gravitational heating 


(e.g. iRees fc Ostrikeilfl977l : Khochfar fc Ostrikerll2008l l have 
been considered. 


Galaxy formation is a complex process which involves 
both stars and gas, and therefore, galaxy formation models 
are required to reproduce the observed properties of not only 
stellar but also the cold-gas components of galaxies. The 
observed mass-metallicity relation of galaxies implies that a 
galaxy is not a closed-box system, but evolves with the con- 
tinual cycling of bary ons between galaxies and the IGM (e.g. 
iTremonti et al.ll2004l l. Such a picture of galaxy evolution 
is sup ported by recent c osmological hydrodynamic simula¬ 
tions ((Dave et al. I l2011a| ia. In such a galaxy-IGM ecosys¬ 
tem context, galaxies are described a s a slowly evolving 
system (i.e. ‘qua si-equilibrium’ model; iBouche et al.l l20ld : 
iDave et al.ll2012tl . where the galaxies acq uire gas from cos¬ 
mic large-scale fil amentary structures le.g. lKeres et al.ll2005l : 
iDekel et al.ll200d l which activate star formation, lose gas by 
strong and ubiquitous outflows which regulate star forma- 
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tion, and again acquire the gas returned from the halo (i.e. 
wind recycling: lOppenheimer et aI.ll2M^ . Therefore, inves¬ 
tigations of the evolution of the gas component in galaxies 
are essential to understand galaxy formation. 

Observations of cold-gas components in higher-redshift 
galaxies are currently underway. The cold-gas components, 
atomic hydrogen gas (Hi) and molecular hydrogen gas 
(H 2 ), are traced by the 21-cm hyperfine structure line 
and the 3-mm line from the carbon monoxide (CO) ro¬ 
tational transition, respectively. Although most Hi emis¬ 
sion studiej^ in high redshift have been limited to the 
z ~ 0.1 — 0.2 re gime due to instrumental frequency coverage 
and sensitivity ( Catinella et al.||2008l: iR eudling et al.|[20i^ : 


iFernandez et al. 20131 : Rhee et al.ll2013[ l. emissions from CO 

have been detected farther out to 2 ~ 6; these emissions 
have b een used to measure t he H 2 mass in the high -2 uni¬ 
verse llCarilli fc Waited l2013l . and references therein). Re¬ 
cent CO surveys of both local and high -2 galaxies have re- 
vealed molecular gas fractions as a function of stellar mass 
llSaintonge et al.ll2011al : [Tacconi et al.|[201^ . In addition to 
these direct measurements of cold gas, newer advances in 
the observational studies of cold-gas components are also 
expected from other indirect methods base d on the star- 
form ation-rate-molecula r-gas relation (e.g. [Popping et al 


2012) , dust continuum (iMagdis et all l2012ain'Santini et al 


2014 ; iscoville et al.l [201^ 2015|), and the spectral features 


in the optical band ( Morokuma-Matsui et alll2015fl . These 
recent developments, as well as observations made with 
the newest generation of radio and sub-mm instruments 
such as ALMA (Atacama Large Millimeter/submillimeter 
Array), SKA (Square Kilometer Array), ASKAP (Aus¬ 
tralian Square Kilometre Array Pathfinder), MeerKAT, and 
NOEMA (Northern Extended Millimeter Array) are ex¬ 
pected to provide us important constraints on galaxy for¬ 
mation theories. 

Recent cosmological galaxy formation models con¬ 
tain sufficient physical ingredients (e.g. gas cooling, 
H 2 formation, star formation, and feedback) such 
that direct comparison with observations is now pos¬ 
sible. Although the statistical properties of galaxies 
have been studied w i th semi-analytical models (SAMs ; 
e.g. i^uffmann et al.| 19991: Sornervill e fc Primacld 1 19991 : 


ICole et al.ll2000l : iNagashima et al.ll2005ll . recent large-scale 
cosmological A^-body/hydrodynamic simulations also al¬ 
low us to statistically study galaxy formation and evolu- 


low us to statistically study galaxy tormation ana evolu¬ 
tion (e.g. Schwejt_^ l2010l. I 2 OI 5 I: lOkamoto et al.l I20l4 


[Thompson et al. 20l4 IVogelsberger et 2014( 1. In addi¬ 
tion, some of these theoretical studies have implemented 
the equilibrium or non-equilibrium formation of H 2 , as 
well as star fo rmation from such H 2 gas, in cosmologi¬ 
cal simulations | Gnedin et al.||2009|: Christensen et aDl2012l : 


lThompsoii_et al.l I20l4 iTomassetti et a l.||2015|) and SAMs 
llLagos et al.l 1^1 Ibl: IFu et al.l 20121 : Popping et al.l I20l4 
ISomerville et al.l 2015lL 

However, these theoretical models still require many 
‘sub-grid’ parameters for star formation and feedback mod¬ 


els due to limited spatial/mass resolutions and a lack of the 
full understanding of related physical processes (for m ore 
details, see a recent review bv ISomerville Davdl2()l4 L As 
described above, there are fewer studies on the evolution 
of the cold-gas component of galaxies than those on stel¬ 
lar components. Recent developments in both observational 
and theoretical studies on the galactic cold-gas components 
provide us a unique opportunity to compare theoretical and 
observational results and to obtain the constraints for theo¬ 
retical models. 

In this study, we investigate as to whether the current 
galaxy formation models reproduce the observed properties 
of the cold-gas components of galaxies by comparing the ob¬ 
served and theoretically predicted evolutions of the molecu¬ 
lar gas mass fraction with respect to the total baryonic mass 
as a function of stellar mass. The outline of the paper is as 
follows. We first investigate the redshift evolution inferred 
from the observational data in section [21 In section O the 
observed redshift evolution is compared with the model pre¬ 
dictions. Further, we discuss feedback models implemented 
in galaxy formation models in section |T] Finally, we sum¬ 
marise the study and comment on possible future prospects 
in this direction in section [5] 


2 OBSERVED EVOLUTION OF MOLECULAR 
GAS FRACTION 


In this section, we investigate the observed molecular gas 
fraction as a function of stellar mass. The molecular gas 
mass fraction is defined as 


fmol 


Mn,oI -h M* ’ 


( 1 ) 


where M* and Mmoi denote the stellar mass and molecular 
gas mass of a galaxy, respectively. We combine the literature 
data of M* and Mmol and calculate the molecular gas mass 
fraction /mol of various galaxies. The literature data of Mmol 
include both direct and indirect estimations. In the following 
sections, we first describe the literature data used in our 
study (sections 12.II and 12.21) and we subsequently investigate 
the redshift evolution of /mol as a function of stellar mass 
(section 12.311 . 


2.1 Measurement of molecular gas mass from CO 

CO data were retrieved from nine studies (|Leroy_et_^|20()8|; 


DadHi_et^lJ[2010yGeachet^lJ[201iySamtoiigeet_alJl201^ 


Bauermeister_et_ajJ^ 20131: iTacconi et al.l 201^ Boselhet_^4 


2014 iBothwell et al.l I20l4 iMorokuma-Matsui et al.l l2oi5n 


to cover a wide range of redshifts and stellar masses. We 
adopted the Galactic conversion factor of aco = 4.35 
M(7)(K km s~^ pc^) ~^ (including the contribution of helium; 
iBolatto et al.|[^13ll for all the sample galaxie^ considered 
in this study to estimate the molecular gas mass as 


Mmol = aco Leo I 


( 2 ) 


^ The absorption line observations of neutral hydrogen are also 


used to test galax 

Y formation theories i 

Rahmati et alj|2013allb 

Rahmati Schave 

I 2 OI 4 !: iRahmati et al. 

I 2 OI 5 I: iBird et alj 12014 

Berrv et al.ll2014'). 


^ The starburst galaxies, which require different acOi are ex¬ 
cluded in this study even though they were included in the origi¬ 
nal studies. 




































































































Redshift evolution of M* * — /moi relation m 0 < z < 2 3 


Name 

Redshift z 

Af* (IOIOMq) 

SFR (Mq/ yr) 

Symbols in figure m 

Reference 

LerOS 

0.0007 - 0.004 

0.001 - 8.0 

0.7-4.1 

open suquare 

Leroy et al. (2008) 

Bosl4 

0.0035 - 0.006 

0.03 - 13.1 

0.01 - 6.0 

cross 

Boselli et al. (2014) 

Botl4 

0.01 - 0.03 

0.03 - 1.0 

0.16 - 4.0 

open pentagon 

Bothwell et al. (2014) 

Saill 

0.025 - 0.05 

1.0 - 32 

0.07 - 40 

dot 

Saintonge et al. (2011a) 

Baul3 

0.05 - 0.3 

4.0 - 30 

3.4 - 88 

filled suquare 

Bauermeister et al. (2013) 

Mori 5 

0.1 - 0.2 

4.0 - 20 

8.5 - 48 

filled circle 

Morokuma-Matsui et al. (2015) 

Geall 

0.4 

4.1 - 11 

28 - 62 

open circle 

Geach et al. (2011) 

DadlO 

1.5 

3.3 - 11 

62 - 400 

open triangle 

Daddi et al. (2010) 

Tad 3 

1.0- 1.5, 2.0-2.3 

0.6 - 17 

26 - 480 

filled triangle 

Tacconi et al. (2013) 

Popl2 

0.5 - 2.0 

0.01 - 100 

- 

dash-dot lines 

PoDDing et al. (2012) 


Table 1. Summary of samples from the literature. Estimations of M* and SFR in each study are described in section [2. II 


where L'qq represents the luminosity of the ^^C O(J =1 — 0) 
line in units of K km s ~ ^ pc^ . The CO data o flLeroy e t al.l 
ll2008li. Bothwel l et al.l (120141 '). iTacconi et al.l ll2013i) . and 
iDaddi et al.l ( Mid i are not based on ^^CO(J = 1 — 0) obser¬ 


vations but based on ^^CO(J = 3 — 2) or ^^CO(J = 2 — 1). 
The conversions to *^^00(7 =1 — 0) intensity assumed in 
these studies are described in the following lists. We discuss 
effects of variable aco on /mol in section 13.21 The literature 
data of the CO observations are summarised in Table [T] We 
also describe the stellar mass and star formation rate (SFR) 
estimations of these literature data in the following list. 


• iLerov et al.l ll2008tl : Leroy et al. combined the literature 

data of the CO maps of local galaxies (3 < D < 30 Mpc) 
including spiral an d dwarf galaxies, ^^CO(J = 2—1) data 
from HERACLES llLerov et al.ll2008ll and ^ ^CO(J = 1-0) 
data from BIMA SONG (iHelfer et al.lliooil . They assumed 
a ^^CO(J = 2 — 1) to ^^CO(J =1 — 0) intensity ratio of 
0.8. The stellar mass and SFR of the sample galaxies range 
over 10^'^ — 10^°'® Mq and 0.7 — 4.1 Mq yr“^, respectively 
and thw are the so-called main sequence of star-forming 
galaxieqj. Leroy et al. derived stellar mass using Spitzer’s 
3.6-^m data and SFR using GALEX FUV and Spitzer’s 24- 
/rm data. _ _ 

• iBoselli et al.l (l2014l l: Boselli et al. observed ^®CO(J = 
1 — 0) emission data of 59 late-type galaxies of the Herschel 
Reference Survey, which is a survey of the complete A'—band 
selected, volume-limited (15 < D < 25 Mpc) galaxies span¬ 
ning a wide range in morphological type and luminosity, 
using the NRAO Kitt Peak 12-m telescope. They combined 
the literature data of CO and provided a CO catalogue of 
225 out of the 322 galaxies of their complete sample with 
stellar mass of 10®'® — 10^^ ® Mq and SFR of 0.01 — 6 Mq 
yr ~^. The stellar mas ses of their samples were estimated 
in ICortese et ahl ll2012ll from f—band luminosities using the 
g — i colo ur-dependent s tellar mass-to-luminosity ratio rela¬ 
tion from IZibetti et al.l (l2009l l . The sa mple galaxies are the 
main sequence of star-forming galaxies (ICortese et al.ll2014l : 
ICiesla et al. 20ML 


Bothwell et al.l (l2014h : Bothwell et al. presented the 


hrst data release of a ^®CO(J = 2—1) survey with the 12-m 
telescope of Atacama Pathfinder Experiment (APEX) for 


® It is known that st ar-forming galaxies f orm a distinct sequence 
of SFR with M*(e.g. iNoeske et al.ll200^ . The galaxies on this 
sequence are called ‘main sequence’ galaxies. 


nearby dwarf galaxies at 0.01 < 2 < 0.03 with the stellar 
mass and SFR ranges of 10®'® — 10^® Mq and 0.1 — 4.0 Mq 
yr“^, respectively. They limited to the samples with metal- 
licity of 12-1-log (O/H) > 8.5 for aco to be similar to or lower 
than the Milky Way value. They assumed a ^®CO( J = 2—1) 
to ^®CO( J =1 — 0) intensity ratio of unity to estimate L'qq. 
The stellar masses and SFRs of their sample galaxies were 
provided by the Max Planck Institute for Astrophysics-John 
Hopkins University (MPA-JHU) group. The stellar masses 
of the galaxies were derived by htting the SDSS ugriz pho¬ 
tometry to models spanning a wide range of star formation 
histories. 


• ISaintonge et al.1 ll2011ah : Saintonge et al. conducted a 
CO survey of local galaxies with stellar masses of (1 — 30) x 
10^® Mq and SFRs of 0.07 — 40 Mq yr“^ with the 30-m 
telescope at the IRAM facility (CO Legacy Data base for 
the GASS survey, COLD GASS), and their targets were se¬ 
lected such that the resulting stellar mass distribution was 
roughly flat. In their study, the stellar mass and SFR of 
galaxies were estimated by fitting the GALEX (NUV and 
FUV) and SDSS photometries withMhejBnrxunl^Charl^ 
(1200311 population synthesis code (ISaintonge et al. 2011bll . 


Their sample consists of the star forming galaxies and pas¬ 
sive galaxies. 

• iBauermeister et al.l (l2013ll : Bauermeister et al. targeted 
31 star-forming galaxies with stellar masses in the range of 
(4 — 30) X 10^® Mq and SFRs in the range of 3 —90 Mq yr~^ 
at z ~ 0.05 — 0.5 and detected CO emissions from 24 galaxies 
using the Combined Array for Research in Millimeter-wave 
Astronomy (CARMA) facility (Evolution of molecular Gas 
in Normal Galaxies, EGNoG). The sample galaxies in the 
ranges of z ~ 0.05 — 0.32 from whic h GO emission was de - 
tected were drawn from SDSS DR7 (lAbazaiian et al]l2009l l. 
The stellar mass and SFR of their sample galaxies were pro¬ 
vided by the MPA-JHU group. This estimate was found to 
be comparabl e with the one obtaine d from stellar absorp¬ 
tion features (iKauffmann et aI1l2003h . Further, SFRs were 
derived by htting the flux es of no less than five emission lines 
dBrinchmann et ^l2004h . In this study, we only use the data 
of the normal star-forming galaxies of their sample. 

• iMorokuma-Matsui et al.l (l2015fl : Morokuma-Matsui et 
al. conducted GO observations of 12 galaxies with stellar 
masses of (4 — 20) x 10^^ Mq and SFRs of 8.5 — 48 Mq yr“^ 
in the z ~ 0.1 — 0.2 regime using the 45-m telescope at the 
Nobeyama Radio Observatory, and they det ected CO emis- 
sion from 8 galaxies. Similar to the case of iBothwell et all 
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||2014i and iBauermeister et al.l (l2013ll . the stellar masses 
and SFRs of their samples were also obtain ed from the MPA - 
JHU group and drawn from SDSS DRIO (|Ahn et al.ll2014h . 
They showed that their samples correspond to star-forming 
galaxies at the redshift. 


• iGeach et al.l ([20l3): Geach et al. observed seven galax¬ 
ies selected based on the basis of the 24-/rm emission galaxies 
with masses of (4 — 11) x 10^^ Mq and SFRs of 28 — 62 Mq 
yr“^ at z ^ 0.4 using the IRAM Plateau de Bure Interfer¬ 
ometer (PdBI) and detected CO emission from hve galaxies. 
The stellar masses of the ir samples were est imated by htting 
BVRIJK photometry (|Moran et ahllioOTh to model spec¬ 
tral energy distribu tions (SEDs) using the KCQRRECT soft¬ 
ware package, v4. 2 l|Blanton fc Roweigl2nn7l ~l. The SFRs of 
their samples wer e derived from the far-infrared luminosity 
dKennicu td 1998 3) est imated from the 7.7-/rm-line luminos¬ 
ity ( Geach et al.ll200^ . 

• IPaddi et al.l (l201Cll l: Their sample galaxies comprised 
the star-forming galaxies with a stellar mass range of (3 — 
11) X 10^° Mq and SFR range of 62 — 400 Mq yr“^ in the 
1.4 < 2 < 2.5 regime class ified by using the BzK colour 
criterion (IPaddi et al.ll20()4 1 with Spitzer imaging detection 
at 24 fj,m. They observed CO(J = 2 — 1) of hve galax¬ 
ies with IRAM PdBI and estimated Lqq by assuming a 
CO(J = 2 — l)-to-CQ(J = 1 — 0) conversion factor of 0.86 
l|Pannerbauer et al.ll200!ih . The stellar mass was estimated 
via an empirical method with BzK photometry, which was 
calibrated with the mass estimates from m ulticolor pho¬ 
tometry and SEP htting llPaddi et al.ll2003l . The SFRs of 
their samples were the averages of three SFRs estimated 
from dust- corrected UV, 24- fim emission, and 1.4-GHz ra¬ 
dio huxes llPaddi et al.ll20fol l. 

• iTacconi et al.l (l2013ll : Tacconi et al. conducted the 
largest CO survey of 2 ~ 1 — 2 galaxies with stellar mass of 
(0.6 — 17) X 10^° Mq and SFR of 26 — 480 Mq yr“^ with 
the IRAM PdBI and detected CO emissions from 52 galax¬ 
ies (IRAM Plateau de Bure HIgh-z Blue Sequence Survey, 
PHIBSS). The stellar mass was estimated from SEP htting. 
The SFR estimate was based on the sum of the observed UV- 
and IR-luminosities for the 2 ~ 1.2 samples lIWuvts et al.l 
l201ll 'l. and the ex tinction-corrected Hq luminosities for the 

2.2 samples jKennicutt jj98a|;_ TOrster Schreiber et al.l 


I 2 OO 9 I : iMancini et al.ll201ll : IWuvts et al.ll201ll l. Their sample 
follows the main sequence of star-forming galaxies at the red- 
shift. They observed the luminosities of the ^^CO(J = 3 — 2) 
lines of galaxies and converted these data sets into the 
luminosity of ^^CO{J = 1 — 0) lines using a ^^CO(J = 
1 - 0)/^^CO(J = 3-2) ratio of 2. 


All the studies except for iLerov et al.l (|20 08ll assum ed the 
Chabrier initial mass function (IMF. [Chabrieil[2003ll for es¬ 
timati ng the stellar mass and SFR of ga laxies. Lerov et al.l 
ll2008tl used Kroupa IMF (lKroupall200ll ~). We did not cor¬ 
rect for this IMF diherence, since the difference in the 
stellar mass and SFR estimations based on Chabrier and 


Krou pa IMFs is small compar e d to the other unc e rtain¬ 
ties (Ichomiuk fc Povichl 1201 ll : iKennicutt fc Evans! l2012l : 
iMadau fc Pickinson 12014 1. 


2.2 Estimation of molecular gas mass from SFR 


We also use the molecular gas mass estimated indirectly 
from optical studies, since the CO measurement is limited 
to galaxies with a relatively narrow range of stellar mass 
(-^ 10^° — 10^^ Mq ) j part icularly at high redshift. 

IPopping et al.l (120121') estimated the .fmoi values of the 
COS MOS (Cosmological Evolution Survey; IScoville et al.l 
l2007h galaxies using the molecular gas mass values inferred 
from the Kennicutt-Schmi dt (K-S) law and the pressure - 
based H 2 formation model (iBlitz fc Rosolowskvll2004l . l2006l l. 
They presented the fitting formulae of the M* — /mol rela¬ 
tions for star-forming galaxies and all galaxies brighter than 
Iab ~ 24 mag. For the star-forming galaxies, the fitting 
formula is given by 

= exp(logM*-A)/R + l’ 


where 

^ = 6 . 15(1 + ^] ’ -B = 1.47(1 + 2)-^•^^ (4) 

and that for the whole galaxies brighter than 24 mag in Iab 
is 


/mol,P12 


1 

exp(logM* -C)/P-f 1’ 


(5) 


where 

C = 6.33(l + ^J , P = 1.71(1+ 2)-i-^^ (6) 

Here, we mention that IPopping et al.l (l2012l l showed 
that their Iab < 24 mag sample is complete (100 %) for the 
range of 10®-^ Mq < M* < 10“-’' Mq at 0.5 < 2 < 0.75 with 
respe ct to the SMF derived in the same field dllbert et al.l 
l2010lf . However, it is almost evenly < 50% for all stellar 
mass ranges at 2 > 1.25. In addition , if we assume the fit¬ 
ting functions of IPopping et al.l (l2012f) as accurate, the /mol 
values of galaxies with stellar masses > 10^^ Mg becomes 
larger at lower redshifts, but the statistics with regards to 
higher -2 massive galaxies in the htting is poor. Thus, it is 
required to closely consider these observational limitations 
when we compare observations with theoretical studies. 


2.3 Stellar-mass-dependent evolution of molecular 
gas mass fraction of galaxies 

In this section, we investigate the stellar mass dependence of 
/mol evolution using the observational data. Figure [T] shows 
the redshift evolution of the fraction /mol as a function of 
stellar mass. Although the stellar mass range of galaxies 
with CO measurements is narrow at high redshift (2 > 0.1), 
we can observe the stellar-mass dependence where the mas¬ 
sive galaxies have smaller /moi values than their less mas¬ 
sive counterparts at 2 ~ 0 — 2. In all redshift ranges, the 
observed M* — /moi.co relations are well reproduced by the 
Al* —/moi,pi2,SF relations in which the molecular gas mass is 
estimated from the indirect method. IPopping et al.l (l2012f) 
have already reported that the M*-/moi,pi 2 ,SP relation suc¬ 
ceeds in repro ducing the observed M *-/moi.co relations of 
local galaxies dSainto nge et al.ll2011ah and 2 ~ 1 — 2 sam¬ 
ples (ITacconi et al.ll2010ll . Here. we found that their M*- 
/moi,pi 2 ,SF relation also reproduces the observed M*-/moi,co 
relations of galaxies at 2 ~ 0.1 — 1. 
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(a) 2<0.1 (b) 0.1<2<0.5 (c) 1.0<2<1.5 (d) 2.0<z 



Bosl4 

S. 

Ler08 
B.iul3 
Geall 
Moris 
DadlO 
Tael ^ 
Popl2 
Popl2-SF 


M* [M.] M* [M. ] M* [M.] M* [M. ] 


Figure 1. Molecular gas mass fraction /mol a function of stellar mass and redshift (a: ^ < 0.1, b: 0.1 < z < 0.5, c: 1.0 < 2: < 1. 5 
and d: 2.0 < z). The symb o ls indi cate CO data fro m the literature ('cross : iBoselli et alj (120141). open-square: iLeroy et al.l 1I2OO8I) . 
small-circle: ISa intonge et al.l ' 2011al) , open pentagon: iBothwell et~a 0 ESI), filled-square: i Ba^rmeister I] 1I2OI3I). filled-circle: 
lMomkuma-Matsu^tal.1 (12015 '). open-circle: ICeach et a iJ iaoiil) , open-triangle: iDaddi et al.l J2OI0I') . and filled-triangle: iTacconi et alj 
(|2Q13l ')). The grey and black dot-dash lines indicate M* — /mol relations based on the K-S l aw and pressur e -base d H2 formation models 
for the star-forming galaxies (/mol,Pl2,SF) ^^d whole samples (/mol,Pl2); respectively in [Popping et al] (120121) . The redshifts of the 
equations for /mol,Pi2,SF /moi,Pi2 0.025 for (a), 0.3 for (b), 1.25 for (c), and 2.2 for (d). The relationship between the symbols 
and references is summarised in Tabled The stellar-mass-dependent /mol evolution, where massive galaxies tend to reduce /mol iii earlier 
epochs compared to their less massive counterparts, can be clearly observed (see section [23] for more details). 


There is a clear stellar-mass dependence in the /moi evo¬ 
lution of star-forming galaxies where less massive galaxies 
tend to show a decr e ase in /moi over time as reported in 
IPopping et al.l l|2012l . l2015l ). Here, we remark that recent 
dust-based studies have also suggested that the stellar-mass 
dependence of the cold gas mass fraction (/gas = MgTs+ m. ’ 
where Mgas denotes the total cold gas mass including the 
atomic and molecul ar gas) shows a sim ilar stellar-mass de¬ 
pendent evolution (ISantini et all 1^141 ). Therefore, studies 
based on two independent methods (CO or dust) both sug¬ 
gest the stellar-mass dependent evolution of the cold gas 
fraction of galaxies. 

We examine the stellar-mass dependent evolution of 
/mol quantitatively. We must consider the stellar mass 
growth when discussing the /mol evolution. If we adopt the 
stellar mass g rowth of star-forming galaxies as reported in 
iLeitnerl (12012 ). the redshift evolution of the M* — /moi,pi2,SF 
relation suggests that galaxies with stellar masses of 10^° 
Mq, ® Mq, and 10^^ M© at 2 = 0 exhibit decreases by 
factors of 8 (from 80 % to 10 %), 6 (37 % to 7.0 %), and 4 
(26 % to 6.1 %) in their molecular gas mass fraction from 
z = 1 to z = 0, respect iveljQ. Similarly, if we use the stellar 
mass growth of galaxies based on the abundance matching 
met hod (i.e. not biased with respect to star-forming galax¬ 
ies; [Con^^T^^chsie^Iiool), the M* — /moi,pi 2 evolution 
suggests that the galaxies with stellar masses of 10^'^, 10^'^'®, 
and 10^^ Mq at 2 = 0 exhibit molecular gas-fraction de¬ 
creases by factors of 14 (from 96 % to 6.8 %), 8 (14 % to 
4.2 %), and 1 (4.9 % to 3.6 %) from z = 1 to z = 0, respec¬ 
tively. Therefore, the /moi evolution of less massive galaxies 
in the range of 2 = 0 — 2 is more drastic than that of massive 
galaxies, and the massive galaxies appear to have depleted 

* * He re we use the M* — /mol,Pl 2 ,SF relation because iLeitnej 
1 I 2 OI 2 I ) estimated the growth rate of stellar mass by assuming that 
present-day star-forming galaxies were on the ‘main sequence’ of 
the star-forming galaxies. 


relatively large fract ions of their molecular gas by 2 : = 10 
dPopping et ahllioisl) . 


3 COMPARISON WITH THEORETICAL 

STUDIES 

In this section, we compare the observed M* — /moi relation 
with predictions of theoretical studies constructed to repro¬ 
duce the observed SMF. The theoretical studies we used in 
this study are briefly summarised in section [3Tl We demon¬ 
strate a discrepancy in the M* — /mol relation between the 
observations and theoretical predictions in section 13.21 

3.1 Galaxy formation models in the literatures 

Recent galaxy formation models take into account star for- 
mati on from molecular gas, as suggest ed by observations 
fe.g. iBigiel et al.ll2008l : iLerov et ^2008h . These theoretical 
models need to consider both star formation and molecule 
formation. Since star formation is an unknown process and 
the current cosmological numerical simulations of galaxy for¬ 
mation cannot resolve the spatial scale of molecular clouds, 
the galaxy formation models treat star formation as a ‘sub¬ 
grid’ model (see below). 

For H 2 formation, most hydrodynamic simulations and 
SAMs adopt simplified methods such as ‘pressure-based’ 
and ‘metallicity-based’ models, although some cosmolog¬ 
ical hydrodynamic simulations solve the non-equil ibrium 
chemical reactions of hydrogen (iGnedin et al.l I 2 OO 9 I ). The 

® .fm nl evolution farther out to ^ = 3 has been also investi¬ 
gated jMagdis et al.ll2012bl : ISaintonge et al]|2013h . iMagdis et alj 
ll2012d) detected a COfJ = 3 — 2) emission from an infrared lu¬ 
minous Lyman break galaxy at 2 = 3.2 with a stellar mass of 
2 X 10^^ Mq. The estimated /mol is 36 % suggesting that massive 
galaxies have flat evolution of /mol least in 1 < 2 < 3. 
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(a) z<0.1 (b) 0.1<2<0.5 (c) 1.0<z< 1.5 (d) 2.0<z 



Figure 2. H 2 gas mass fraction ^ ^ function of stellar mass and redshift (a: 2 < 0.1, b: 0.1 < 2 < 0.5, c: 1.0 < 2 < 1.5 and d: 
2.0 < 2 ). The black and grey dot-dash lines indicate the M* — relations based on the K- S law and pressure-ba sed H 2 formation 
models for the star forming galaxies (/h 2 ,pi 2 ,Sf) whole samples (/h 2 ,pi 2 )) respectively in iPopping et akl (120121) . The redshifts of 
/h 2 ,pi 2 ,SF Q-nd /h 2 ,pi 2 are 0.025 for (a), 0.3 for (b), 1.25 for (c), and 2.2 for (d). [Upper panel] Comparison with the predictions from 
cosmological A ^-body/hydrodynam ic simulations (SIMs; dark-gree n, dark-grey, and b lue dotted lines: high-, middle -, and low- resolutio n 
calculations of iDuffv et al.l ll2012l') . respectively, solid green lines: iGenel et a y ll2014l) . and light-blue dashed lines: [Lagos et ^ ll2015h ). 
ICenel et al.l l|2014^ provided both lower and uppe r limits for for galaxie s at 2 0. [Lower pane l] Comparison with t he predictions from 

semi -analytic model s (SAMs; pink dotted lines: [Obreschkow et aU (120091) . purple dashed lines: [bagos et ^ (l2011al) . light-purple solid 
line: IFu et alJ ll2ni2h . orange dot-dot-dashed line: IPoppin^^ al.l ll2014^ 'l. In the graph legend, ‘P’ denotes the pressure-based H 2 model, 
‘Z’ denotes the metallicity-based H 2 model, and ‘p’ indicates that the H 2 model was implemented as a post-process. The relationship 
between the each symbol and the reference is summarized in Table [2] On each plot, we indicate the theoretical predictions for /hj at 
z = 0 in (a), z = 0.5 in (b), z = 1 in (c), and z = 2 in (d). Theoretical studies ‘qualitatively’ reproduce the stellar mass dependence of 
/hj evolution but there is a ‘non-negligible’ discrepancy between the observations and theoretical studies (see section [3.21 for details). 


Label 

Method 

H 2 prescription 

Star formation model 

Feedbacks 

Line Styles in Fig. [2] 

Reference 

Dufl2 

SIM 

Pressure* 

Schave Dalla Vecchia ('2008) 

SN 

dot 

Duffv et al. (2012) 

Genl4 

SIM 

_* 

Sorineel & Hernauist ('2003! 

SN & AGN 

solid 

Genel et al. f2014) 

Lagl5 

SIM 

Metallicityt * 

Schave &; Dalla Vecchia ('2008) 

SN & AGN 

dash 

Laeos et al. (2015) 

Obr09 

SAM 

Pressure* 

Kennicutt (1989) 

SN & AGN 

dot 

Obreschkow et al. ('2009) 

Lagll 

SAM 

Pressure 

Bigiel et al. 12008) 

SN & AGN 

dash 

Laecos et al. f2011a.b) 

Ful2 

SAM 

Metallicity 

Bigiel et al. 12008) 

SN & AGN 

solid 

Fu et al. 120121 

Popl4 

SAM 

Metallicity 

Bieiel et al. ("2008) 

SN & AGN 

dot-dot-dash 

PoDDine et al. ('2014) 


Table 2. Summary of the t heore tical studies used in Fig. [2] Details of each study are described in section [3.11 tThe H 2 prescription 
is from lOnedin fc Kravtso^ |2oO), which provided a fitting formula to their galaxy simulations incorporating non-equilibrium chemical 
network of hydrogen and helium and non-equilibr ium cooling and heating rates. The other ‘Metallicity’ models use an analytical formula 
based on the chemical equilibrium model of H 2 llKrumholz et al]l2n09al l. *H 2 formation (partitioning) models were implemented as a 
post-process. 


pressure-based model is an empirical method based on the 
observed relation between the molecule fraction of cold 
gas and the surface densities of stellar and gas compo¬ 
nents in nearby disc g alaxies (iBlitz fc Rosolowskvl |2004| . 
I2OO6I : iLerov et ^l 2008 l ). In contrast, the metallicity-based 
model uses an analytical formulae of the molecule fraction 
as a function of the gas density, m etallicity, and radiation 
held llKrumholz et al.l 20081 . [2009a^ . It has been shown that 
both these types of models suitably reproduce the observed 


molecular fraction of cold gas well llBl itz fc Rosolowsk'^ 
|2004| . I2OO6I : iKrumholz et al.|l2008l : iFumaealli et al.ll201(]|) . 


In this study, we compare the observations to the 
predicted M* — /mol relations from three cosmologi¬ 
cal V-body/hydrodynamic simulations (S IM; iDuffv et al.l 


I 2 OI 2 I: Genel et al. 12014: iLaeos et al. 

2015ll and four SAMs 

llObreschkow et al.l 

20091: iLaeos et al. 

l2011al:lFu et al.ll2012l: 

IPoDDinsr et al.ll2014 

) of galaxy formation. We summarise the 

theoretical studies used in this study 

in Table [2] 
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• iDuffv et alJ ll2012l . SIM): Duffy et al. implemented the 
press ur e-based H 2 forma tion model (iBlitz fc Rosolowskvl 
l200fil: iLerov et all l2008h in the calculations of the 
OverWhelminRly L arge Simulations (OWLS) project 
dSchave et al.lHoiOl 'l as a post-process. They determined 
the free parameters of the H 2 formation model according 
to recent observations dLerov et al.l 1200^ 1. Their estimated 
H 2 mass included the contribution of helium. The simu¬ 
lations in th e OWLS include radi ative cooling (10^ K < 
T < 10® K IWiersma et all l2009l ~). photoionisation heat- 
ing from cosmic UV backgrou nd radiation, star formation 


SNe fi.e. bv ‘kick’: Dalla Vecchia & Schavel 

2008 

) and ther- 

mal feedback from AGN (|Booth & Schavel 

2009 

). Star for- 


ll2008tl . in which they analytically recast the K-S law as a 
function of pressure rather than density, while assuming a 
self-gravitating disc. Duffy et al. provided results for the dif¬ 
ferent co-moving box sizes: 25 h~^ Mpc (high resolution), 50 
h~^ Mpc (middle resolution), and 100 h~^ Mpc (low resolu¬ 
tion). The simulations with the smallest box size were run 
to 2 = 2 and t hose with the large r box sizes were run to 
2 = 0. Although iDuffv et al.l (l2012l ~) examined both cases of 
with and without AGN feedback, we adopted the former case 
since Duffy et al. already reported that the latter results in 
a la rge discrepancy bet ween the simulation and observation. 

• iGenel et al.l (l2014l . SIM) Genel et al. did not implement 
any H 2 formation models but assumed that one-third of the 
neutral gas of the galaxy is in the molecular phase to derive 
the lower limit f or /h., of 2 ~ 0 galaxies f ormed in the Illus- 
tris simulation dVogelsberger et al.ll2013 L Their estimated 
H 2 mass does not include the contribution of helium. They 
also assumed that all the neutral gas in a galactic disc is in 
the molecular phase to derive /h 2 of high -2 galaxies and the 
upper limit of 2 ~ 0 galaxies. The Illustris uses a moving 
mesh code, AREPO llSpringelll2010l i and includes radiative 
cooling, heating from time-dependent UV backgr ound radi¬ 
ation, star formation dSpringel fc Hernauistll^03l L mechan¬ 
ical feedback from SNe and thermal/mechan i cal/radiative 
feedbacks from AGNs ( ISpringe l et al.l l2005al : ISiiacki et al.l 
I2OO7I : IVogelsberger et al.ll2013il . 

• iLagos et al. ll2015l . SIM) Lagos et al. implemented 
H 2 formation models as a post-process in EAGLE (Evo¬ 
lution and Assemb l y of GaLaxies and their Environ¬ 
ments; ISchave et al.l uOlsIl . which is a cosmological nu¬ 
merical simu lation that u s es a n Wbody/SPH code, 
GADGET-3 dSpringel et al.l I2OO8I L Here, H 2 formation 
is modelled according to two prescriptions depending 
on the local dust-to-gas ratio and t h e interstel l ar ra - 
diation field liGnedin fc Kravtsov! 1201 ll : lKrumhoi3 I 2 OI 3 II . 
Their estimated H 2 mass does not include the con- 
tribution of helium. They considered r adiative cooling 
dWiersma et al.l I2OO9I : IScliave et al.l l2015h . photo-heating 
by cosmic microw ave background, UV an d X-ray back- 
groun d radiation (|llajrrdt_j^^ladau| 1200 il l, star forma¬ 
tion dSchave fc Dalla VecchiaT 2008ll . and the rmal feed¬ 
backs from SNe jPalla Vecchia fc Schavd I2OI2I I and AGN 
llBooth fc Schavell2009l l. We employed their A L — /moi re- 
lation (‘Ref-L100N1504’ model, see Fig. 5 of iLagos et al.l 
2015|) calculated using the H 2 formation model of 


the pressure-based H 2 f o rmation (partitio ning) model 
dBIitz fc Rosolowskvl I2OO6I : iLerov et al.l I2OO8I I as a post- 
process in the s emi-analytical model, L-GALAXIES 
dCroton et al.ll200^ : IPe Lucia fc Blaiz^l2007l l , in which the 
dark matter halo and sub-halo merg er trees are constructe d 
from the Millennium Simulation dSpringel et al.l l2005bll . 
With the L-GALAXIES, they implemented radiative coolin g 
with the cooling function of SuHierlmd&Do£ita|_|l99^L 
phot oionisation heating_|GnedM 2000|: Kravtsov^£_^ 
l2004h . star formation iKa^^nang 

and the feedbacks from SNe and AGN ( Groton et al.l 
2OO6II. We ob t ain t heir A/hj/AI* values from Fig. 9 of 


Popping et al.l d 2015ll and recalculate /h 2 • Here, we mention 
that IPopping et al.l 1 2 OI 5 II did not include the contribution 


from helium. 

• ILagos et al. 1 d2011blfl SA M): Lagos et al. imple - 
mented both pressure -based ([Blitz fc Rosolo wskvl I2OO6I I 
and metallicity-based dKrumholz et al.l 2009 al i H 2 forma- 
tion models in t he semi-analytical model, GALFORM 
(Cole^jt^ I2OOC1I : I Benson et al.l l 2003 l : iBaugh et al.l l 2005 l : 
I Bower et al.l 2006 f l. in which the merger trees of dark mat- 
ter halo were c onstr ucted using the Monte Garlo method 
dLacev fc Cold Il 993 tl . The y considered the shock h eat- 
ing and radiative cooling ( Sutherland fc Dopital 7 99311 in¬ 


side t h e dark ma t ter ha l os, star formation (Lagos_^^_^ 


Bigiel et al 


2011 b: lKennicut '3 Il 998 b l: ILerov et al.l I2OO8I : 

20081 : Krumholz et al. l2009aL and the feedbacks from SNe 


and AGN ( Baugh et al.ll2005 : (Bower et al.ll2006l L We obtain 


their values (pr essure-based H 2 m odel and star 

formation model bas ed on iBigiel et al.|[2008l 'l from Fig. 9 of 
IPopping et al.l (I 2 OI 5 II and recalculate /h 2 . Here, we mention 
that IPopping et al.l ( 2015li did not include the contribution 
from helium. 

• IFu et all (l201Cll . 20lj, SAM): Fu et al. implemented both 
pressur e-based j BIitz fc Rosolo wskvl 1200^ 1 and metallicity- 
based (iKrumholz et akP 2009a|) H 2 f o rmation models in 
the L-GALAXIES ( Croton et al. 200^ De Lucia fc Blaizo^ 
2007|) as is the case with Obreschkow et al. ( 20091 '! . Unlike in 


Obreschkow et ahl (|2009| L H 2 formation in this study is cal- 


Gnedin fc KravtsovI ( 20 111) . 


culated in a self-consistent manner. The mass contr ibution 
from helium was tak en into account (a factor of 1.4. lArnettI 
1 19961: [rT et al. IH). Four kinds of star formation models 
were implemented: 1) ‘B igiel’ (i Bigiel et al.|l200 8l: Lerov et al.1 
I 2 OO 8 ), 2) ‘Kennicutt’ dKennicuttl 19891 . fl )9^ ), 3) ‘Gen- 
zel’ ( Genzel et al.llioiol l and 4) ‘Krumholz’ ( Krumholz et al.1 
l2009bl 'l models. They concluded that the ‘Bigiel’ model re¬ 
sults in better agreement with the observed mass-metallicity 
relation among the four models. We employ their Af* — /mol 
relation calculated using metallicity-based H 2 formation and 
the ‘Bigiel’ star formation model for a comparison purpose 
in this study. 

• IPopping et al.1 (I20l4 SAM): Popping e t al. imple- 
mented pressure-based dBlitz_&_^jOSolows^ I 2 OO 6 II and 
the metallicity-based (iGnedin et al.l l2009ll H 2 forma- 
tion models in the S a nta C r uz semi-analytica l mode l 
dSornerville fc PrimackI 1 19991 : ISomerville et al.l l2008al . 
l2012ll . in which the dark matter halo merger tree was 
constructe d using the Exten d ed P res s-Schechter for¬ 
malism dSomerville fc Kolatd 1 19991 : ISomerville et ahl 
[2008^L The y considered photo ionization heating 
(Gneduj 120001; |Kraytsoy_et_ ^ l2004|). radiative coolin g 


I Obreschkow et ^ d2009l . SAM): They implemented ( Sutherland fc Dopita 19931 : Somerville fc PrimackI ll999l L 
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star formation from molecular gas dBigiel et al.l l2008ll . 
and feedbacks from SNe and AGN. We obtain their 
values (metalli city-based H 2 m odel and star 
formation model bas ed on iBigiel et al.|[2008l 'l from Fig. 9 of 
IPopping et al.l ll 2015ll and recalculate /h 2 • Here, we mention 
that [Popping et al.l ( 2015ll did not include the contribution 
from helium. 

We use the H 2 gas mass fraction /h 2 = ’ "^hich does 

not include the contribution form helium, instead of /mol for 
a comparison between the observations and the theoretical 
predictions, since most theoretical studies considered here 
derived /hj • For the studies presenting /mol instead of /hj , 
we recal culate Mho f rom M moi using a convers i on fac tor of 
1.36 for iDuffv et al.l (l2012fl and Popping et al.l 1 I 2 OI 2 I . indi¬ 
rect method), and 1.4 for IFu et alj 1 I 2 OI 2 II . Hereafter, the 
/h 2 calculated based on equation m and @ is described as 
/h 2 ,pi 2 ,sf and the one calculated based on equation ((5]) and 
is described as /h 2 ,pi 2 . 


3.2 Discrepancy between observations and 
theoretical studies 


Figure [2] compares the observed and predicted M* — /h 2 
relations from theoretical studies. Here, we adopt /h 2 ,pi 2 
(black dot-dash lines in Fig. [2]) as an observational refer¬ 
ence for comparison with theoretical studies since the the¬ 
oretical studies are not limited to star-forming galaxies. It 
should be noted that the theoretical studies did not apply 
any magnitude limits, which were applied to the observa¬ 
tions. This uncertainty is di scussed below. The predicted 
Mi,— fa-, relations from SIMs (iDuffy et al.l20l3: Genelet_^ 
201^ La.gos et al.ll201^ and SAMs llObreschkow et al.ll2009l:^ 


La gos et al.l 2011bl a: et al. I I 2 OIC 1 I. I 2 OI 2 I : IPopping et al.l 

2 OI 4 II are overlaid on the observed relation shown in Fig. 


The theoretical studies qualitatively reproduce the 
trends of the stellar-mass-dependent evolution of /h 2 where 
the less massive galaxies have higher /h 2 values than their 
massive counterparts; however, there is a ‘non-negligible’ 
gap between fur, from observations and theoretical stud¬ 
ies ([Popping et al.ll^OlSl l. The /h 2 predictions from the the¬ 
oretical studies appear to be consistent with the /h 2 ,pi 2 
values in the stellar mass range of ~ 10^° — 10^^ Mq in 
Fig. [2 k. which is consis tent with the results of a recent 
study dbagos et al.[[2015l 'l. However, the predicted val¬ 
ues by the SAMs of less massive (< 10^'^ Mq) are under- 
estimated with re s pect to the observed values except for 
lObreschkow et ^ (l2009[ l. The discrepancy between the ob¬ 
servations and theoretical predictions is still observed even 
if we only consider nearby galaxies (Figs. [2k, b) in the stel¬ 
lar mass range of 10^'^Mq < M* < 10^^ Mq, wherein the 
obs ervational completeness is 100 % (see section [2^ . 

lObreschkow et ^ ([20091 1 predicted higher /h 2 values 
than the observation values for less massive gala xies, even 
thoug h they used the same code (L-GALA XIES) as iFu et 111] 
([ 201 ^ . However. lObreschkow et al.l ( 200^)_Jmplemented H 2 
formation as a post-process whereas Fu et al.l ([2012tl calcu¬ 
lated H 2 formation in a self-consistent manner. This differ¬ 
ence m ay explain the discrepan cy between the two mod¬ 
els, but lObreschkow et ^ (120091 1 claimed that their galax¬ 
ies with a stellar mass of < 4 x 10® Mq were typically 


found in a halo with less than 100 particles, whose merg- 
ing history could be follo wed over only a few time steps 
(iDe Lucia fc Blaizot[[2007l l. 


There is a difference in the mass dependence of /h 2 5 
even though the SIMs and SAMs predict similar red- 
shift evolution of the SMF of galaxies (see Fig. 4 of 
[Somerville fc Davd[20lJ l. At all the redshift ranges (0 < 
z < 2), the SIMs predict more flat M* — /h 2 relations 
than those of the observations (upper panels of Fig. [2]), 
whereas the SAMs predict kinks at the characteristic masses 
of 10^® — 10^^ Mq (lower panels of Fig.[2|). The difference in 
the prescriptions of H 2 formation is not likely the reason for 
this difference between the SAMs and SIMs since there is no 
clear difference among the SAMs depending on whether or 
not H 2 formation is implemented as a post-process (lower 
panels of Fig. mi o r depending on the models (pressure- 
or metallicitv-based: IFu et al][2012k . Thus, the difference in 
modelling star formation and feedback may cause this dif¬ 
ference. We discuss this gap between the SIMs and SAMs in 
section 


In the following paragraphs, we discuss possible causes 
for the /h 2 discrepancy between the observations and the¬ 
oretical studies from an observational point of view. The 
observational studies essentially tend to be biased towards 
bright objects such as gas-rich objects and/or optically 
bright objects if the targets are selected based on large opti¬ 
cal surveys. It has been claimed that the models reproduce 
the observational results relatively well if they take into ac¬ 
count t he sample selection criteria and observationa l sensi¬ 
tivities dKauffmann et aT1[2012l : [Popping et al.[[^014h . These 
studies showed that the predicted /h 2 value approaches the 
observed value if these observational conditions are consid¬ 
ered, but there still seem to be discrepancies between obser¬ 
vatio ns and model predictions (see Fig. 11 in [Popping et ah] 
[2014h . In addition, a gap between a semi-empirical model, as 
a proxy fo r observations, and th eoretical predictions was also 
reported ([Popping et al.[r2015h . This semi-empirical model 
couples a halo abundance matching model with a pressure- 
based H 2 formation model and is not affected by the ob¬ 
servational limitations. However, we must pay attention to 
this observational limitation when comparing the observa¬ 
tions and theoretical studies, especially for the less massive 
and higher-z galaxy samples, which are expected to be more 
biased to bright (gas-rich) objects compared to massive and 
I0W-2: galaxy samples. 


In addition to the observational bias, the uncertainty 
in the CO-to-H 2 conversion factor, aco, also contributes 
an add itional error in the /h 2 estimate. [Narayanan et af] 
(|2012b[ ~l formulated Pco as a function of the gas- 
phase metallicity and l uminosity-weighted CO intensity. 
[Narayanan et all ([2012a[ l showed that the observed /moi 
value with Narayanan’s aco exhibits better agreement with 
the predicted values obtained with cosmological simulations. 
Even though the dispersion is reduced, certain systematic 
stellar-mass-dependent residuals still remain (see Fig. 3 in 
[Narayanan et al.l2012al '). Therefore, the observational causes 
alone may not be able to explain the discrepancy. 
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4 IMPLICATIONS FOR FEEDBACK MODELS 

As described in section there is a gap between the ob¬ 
served M* — /h 2 relations and theoretically predicted rela¬ 
tions. Although some fraction of this gap can be attributed 
to observational uncertainties such as the selection bias and 
variable Qco, the gap still remains. This suggests that the 
gap might originate in the uncertainties of the theoretical 
models. In this section, we focus on the uncertainties of the 
feedback processes in the theoretical modelijf]. The free pa¬ 
rameters for the feedback processes are fixed to reproduce 
the observed SMF even though the SMF is not a direct out¬ 
come of the feedbacks alone, but a consequence of the sev¬ 
eral processes involved in galaxy evolution. This is because 
there is no definitive observational constraint for each feed¬ 
back process. There might also be room for improvement in 
the H 2 and star formation models in terms of the molecule 
fraction and the K-S law, but they are set to meet the direct 
observational constraints for these processes. 

To investigate the stellar mass and redshift dependence 
of each feedback process considered in theoretical studies, we 
construct a n equil ibrium model (EM) based on the study by 
IPave et al.l (l2012l 'l. In the following subsections, we briefly 
describe the EM fsection l4.ll) . and subsequently, we examine 
the implications of the feedback processes through compar¬ 
isons between the EMs and observations fsection l4.2|) . 


4.1 Equilibrium model for galaxy formation 

The EM model is based on the idea that the simulated galax¬ 
ies evolve in a q uasi-equilibrium state (iBouche et al.ll2010l : 
iDave et al.1120121 1: 

M'm = Mi, Mout, (7) 

where Min, Mout, and M* denote the mass inflow rate, mass 
outflow rate, and SFR, respectively. This model assumes 
that the stochastic variations in Min, arising due to factors 
suc h as mergers, are p rone to return galaxies to equilibrium 
fsee lPave et al.lf 2012 l . for details). 

4 . 1.1 Stellar mass and gas fraction 

In this section, we briefly introduce the formula tion of stel¬ 
lar ma ss and gas fraction in the EM according to lPave et al.l 
l| 2012 tl . Integrating the equation M* = Min — Mont, we ob¬ 
tain the stellar mass M*(Mh, z) as a function of the redshift 
z and halo mass Mh. The outhow rate is given as 

Mout = ??M*, (8) 

where the mass lo ading factor ?? is described as 77 = 
^o) (jPave et al. I l 2011 al ig). The inhow rate 
into a galaxy Mn is expressed as 


Min — Mgrav “ Mprev + Mrecyc, 

(9) 

where Mgrav, Mprev, and Mrecyc denote the baryonic inflow 
rate into a galaxy’s halo, rate at which material ends up 

in the gaseous halo of the galaxy, and return 
respectively. These are expressed as 

infall rate. 

Mgrav = fbMh, 

( 10 ) 

-^prev ~ (1 

( 11 ) 

a’/T nV 

1 — az 

( 12 ) 


Here, Mh denotes the halo mass growth rate as given by 


Mh = 25.3 j (l+1.65z)\/lIn>(l + ^)3+nA,(13) 

which is derived fro m cosmological A-body simulations 
dFakhouri et al.l[2010ll , ft, represents the cosmic baryon frac¬ 
tion (0.17), fim the matter density parameter (0.3), Ha the 
dark energy density parameter (0.7), ^ the preventive feed¬ 
back parameter (a product of four feedback parameters, 
CphotoCgravCquenchCwind; see the next subsection for 
explanations on each feedback), and az the metalhcity ratio 
of inflow gas and IS^fl qz = (0.3 - 0 .l 2 )(M*/ 10 ^° Mo)°'^® 
dPave et al.l 12012 1. Substituting these equations into equa¬ 
tion m yields 

Min = -Mgrnv. (14) 

1 — az 

If we substitute equations and d into equation 0 , 
the SFR can be written as 


M* 


CMgrnv 

{1 + r]){l - az)' 


(15) 


In the EM, if M* and M* are known via the procedure 
described above, the gas fraction of a galaxy is given by 


/gas — 


Mgas 

Mgaa -f M* 


1 

l-h(tdepM*/M*)-l’ 


(16) 


where tdep = Mgas/M* denotes the depletion time of the gas 
component by star formation. We use the Hubble time tn 
to parameterise fden = 0.4tH(M*/10^° Mq)“°'® as discussed 
in IPave et al.l d2nid 'l. 


4 . 1.2 Preventive feedback parameters 

Each preventive feedback parametei0 ^ depends on the halo 
mass Mh and redshift z. The EM considers four feedback 
processes as follows: 

• Cphoto, photoionising heating (hereafter PHOTO): An 
ultraviolet background radiation heats the halo gas of galax¬ 
ies. It effectively works in less massive systems since less 


° Recent studies on galaxy formation based on the semi-analytic 
approach have suggested that the gap between the observa¬ 
tions and theoretical predictions decreases if the timescales for 
re-accretion of the ejected gas have a mass- and time-dependence 
jHenriaues et all 1201.31 . l201.^ : iMitchell et al.l l2014l : IWhite et alJ 
I 2 OI 5 II . Even though the physical origin for this treatment is not 
clear, some feedback processes may be attributed to changing the 
re-accretion timescale. 


iDave et al.l ||2012|~) used az = (0.5 - 0.l2)(M*/10l° M0)°'25 
but we use az = (0.3 — to prevent the 

results from being diverged. This difference is because some of 
the equations, which are not fully based on physics o r ab initio 
simulations, were not presented in iDave et al.l l[2o3) (such as 
Mq(z) for equation (1191 1 and Mwind(^) for equation 02011 ): hence, 
we must formulate them by ourselves. 

® In Dave’s model, all feedbacks are modeled in ‘preventive’ man¬ 
ner but some are represented as ‘ejective’ feedback. 
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z=0 z=l z=2 



c 

^hoto 

r 

^quench 

c 

^grav 

r 

^wind 


Figure 3. Preventive feedback parameters as a function of halo mass at 2 = 0, 1.0, and 2.0. The total preventive feedback parameter 
(black solid line) is the product of four preventive feedback parameters, Cphoto dashed lines), Cgrav (gray dotted lines), 

Cquench (gf^'Y short-dotted lines), and CwiND (gray dashed lines). 


massive galaxies cannot confine or acquire heated (pho- 
toionised) gas due to their shall ow potential wells (e.g. 
lGnedinl|20Qd : lOkamoto et al.|[20Q8l ). 


Cphoto = 


^ 1 

3 I Mhalo J 


(17) 


wher e the photosuppression mass, M-y[z) is taken from Fig. 
B1 of lOkamoto et all (|2008l i. 

• Cgrav, virial shock heating during gas accretion (here¬ 
after GRAY): It effectively works in massive systems since 
the sparsity and high virial temperature of the halo gas in 
massive system s prevent the halo gas from accreting onto 
the galaxy ('e.g. iFaucher-Giguere et al1l201ll '). 


Cgrav = MIN 


1,0.47 


l + ^sO.38/ 

4 ) (1012 M 0 


-0.25 


.(18) 


• Cquench, heating up the ISM and/or halo gas and pre¬ 
venting cooling flow from the halo due to supermassive black 
hole (BH) growth (or AGN feedback; he reafter QUENGH) : 
This feedback wor ks in massive systems dCroton et al.l[200a : 
iGabor et al.l[2nnh . 

-1.5 

(19) 


Cquench — 


, 1 / Afhalo \ 

3 Im,(*) j 


where the quenching mass is given by Mq{z) = 
2q12-i-o.3(i-i-z) formulate M qjz) s o as to represent 

the red line in Fig. 1 of iDave et al.l (l2012l l. 

• CwiND, heating of surrounding gas by SNe (hereafter 
WIND): it effectively works in less massive systems since 
ubiquitous SNe wind tends to blow out the gas in less 
massive galaxies due to their shallow potential wells (e.g. 
lOppenheimer et'^l2010l l. 


CwiND 


1 C Mwind(2)\l 

50 V Mhalo )\ ’ 


( 20 ) 


where the wind mass is given by Mwind = 
10i2.5-o.5(i-i-z) ]y|-^ -yyTg formulate Mwind so as to rep¬ 
resent the blue line in Fig. 1 of I Dave et al.l (l2012ll . 


We present the preventive feedback para meters as a funct ion 
of halo mass in Fig. [3] (see also Fig.l of lPave et al.llioi^ . 


4.2 Comparisons between the EMs and 
observations 


We compare the observed M* — /gas relations with the pre¬ 
dictions from the EMs in Fig. 2] Here, we use Popping’s 
Mi, — /gas relation as the observational reference (Popping 
et al. in prep.), which is described below: 


/gas,P12 


1 

exp(logM*-£:)/F-bl’ 


( 21 ) 


where 

E = 9.04(l-f—j , E = 0.53(1-f (22) 

We use the M* — /gas relation because Dave’s EMs do not 
consider H 2 formation. To clarify the contribution from each 
feedback process, five conditions are considered: 1) all feed¬ 
back processes are considered, 2) PHOTO is not considered, 
3) GRAY is not considered, 4) QUENCH is not considered, 
and 5) WIND is not considered. 

From Fig. [I] we can confirm the stellar-mass depen¬ 
dence and observe the degree of influence of each feedback 
process as expected from Fig. [3] For less massive galax¬ 
ies (but > 10®'^ Mq according to observational complete¬ 
ness), the theoretically predicted /gas value approaches the 
observation value at 2 = 0 (Fig. |4^), but it is consider¬ 
ably lower than the observation value even if we exclude the 
PHOTO and WIND feedback processes. On the other hand, 
the QUENCH and GRAY feedback processes work well for 
massive galaxies. The stellar-mass dependence of /gas in the 
high-mass range (but < 10*^^ Mq according to observational 
completeness) appears to be different in that there is a knee 
around M* ~ ® Mq in the EM as opposed to a mono¬ 

tonic decrease in observations; however, it is not conclusive 
because these are consistent within the typical dispersion of 
the observed relation. 

From Fig. 21 we observe that theoretical studies tend to 
predict lower gas fractions than observations, particularly in 
less massive and massive galaxies. Our result suggests that 
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the feedback models implemented in the theoretical stud¬ 
ies are clearly required. However, the physical ingredients 
in these models are still insufficient to correctly reproduce 
the available observations. One requirement is to consider 
a feedback model that quenches star formation but retains 
the cold-gas component in galaxies. 

The EM predicts a stellar-mass dependence of /gas sim¬ 
ilar to that of /h 2 predicted by the SAMs, rather than the 
SIMs; there is a kink at the characteristic stellar mass. The 
Af* — /gas relations predicted by the EMs without (^grav or 
Cqhench become flat, similar to the predictions of the SIMs. 
The observations indicate that the M* — /gas or M* — /hj re¬ 
lations are similar to the predictions of the EM and SAMs at 
higher redshifts, while these relations are similar to the SIMs 
at lower redshifts. The virial shock heating is expected to be 
included in the SIMs naturally, but AGN feedback is imple¬ 
mented as an ad-hoc mod el. If the obser v ed /g as evolution 
in the high-mass range of IPopping et all (120121 ') is correct, 
our comparison suggests a more drastic redshift evolution of 
AGN feedback than the SIMs, i.e. a stronger AGN feedback 
for a high redshift and a weaker one for a low redshift. How¬ 
ever, the details regarding the modelling of AGN feedback 
are beyond the scope of this paper. 

Some recent studies have emphasised the relevance of 
the feedback processes, which suppress star formation while 
simultaneously preserving the cold gas dMakiva et al.l[^014l : 
iMorokuma-Matsui et al.1l20l5ll . Our results showed that the 
EMs and the SAMs tend to underestimate /gas or Jh^ in less 
massive and massive galaxies. There is a larger gap between 
the observations and theoretical predictions for less mas¬ 
sive galaxies in Fig. 0^ (M* — /gas relation) than in Fig. [5^ 
(M* — /h 2 relation). This larger gap in /gas in less massive 
galaxies suggests a lot of ejection of cold gas by SN feedback 
modeled in the theoretical studies. The gap may be reduced 
if we consider the suppression of H 2 formation instead of the 
ejec tion of total cold gas . 

iMakiva et al.l ll20l4 ) implemented a star formation law 
with a feedback depending on galaxy-scale mean dust opac¬ 
ity and metallicity in SAM and succeeded in reproducing the 
faint-end slopes of galaxy luminosity functions at z = 0 with 
a reasonable strength of the SN feedback. In their model, star 
formation in less massive galaxies is suppressed because such 
galaxies tend to have lower metallicity, lower surface density, 
and thus, lower H 2 fraction. For massive galaxies, morpho¬ 
logical quenching is one of the important processes, which 
can que nch star forma t ion ev en in the presence of a gas com¬ 
ponent. iMartig et al.l ll2009li showed that star formation is 
suppressed once the central spheroidal component becomes 
sufficiently massive to stabilise the galactic disc against lo¬ 
cal gravitational instability. Some observational studies sug¬ 
gested that mo rphological quenching is in action in the early- 
type galaxies (iMartig et al.l l2013l : iMorokuma-Matsui et al.1 

I 2 OI 5 I) . 


5 SUMMARY 

Galaxy formation models are constructed so as to repro¬ 
duce the observed SMF of galaxies by blowing out and/or 
heating up the cold-gas component to quench further star 
formation, particularly in less massive and massive systems. 
To investigate whether the current galaxy formation mod¬ 


els also reproduce the observed properties of the cold-gas 
component of galaxies, we compare the observed and the¬ 
oretically predicted evolution of the molecular gas fraction 
with respect to the total baryonic mass as a function of stel¬ 
lar mass (i.e. the M* —/moi or M* —/h 2 relations). Our main 
findings are as follows: 

(i) The evolution of the observed M* —/moi relation shows 
a stellar-mass dependent evolution; massive galaxies have 
largely depleted their molecular gas at 2 = 1, whereas less 
massive galaxies tend to convert molecular gas into stars in 
the regime of 0 < 2 < 1. 

(ii) The current cosmological hydrodynamic simulations 
and semi-analytical models for galaxy formation succeed in 
‘qualitatively’ reproducing the M* — /h 2 relation evolution 
in the regime of 0 < 2 < 2 where less massive galaxies are 
always more gas-rich than massive galaxies. 

(iii) There is a non-negligible gap between the observed 
and theoretically predicted M* — /h 2 relations. The pre¬ 
dicted /mol value is more than two times smaller than the 
observed value at most, particularly in less massive (< 10^° 
Mq) and massive galaxies (> 10^^ Mq). Accordingly, for 
this gap, the use of feedback models could be required, which 
quench star formation but retain the cold-gas component. 

(iv) The observed redshift evolution of the cold gas mass 
fraction in massive galaxies suggests a more drastic evolution 
of AGN feedback than the one implemented in the hydrody¬ 
namic simulations, i.e. a stronger AGN feedback in higher 
redshift (2 > 1) and a weaker one in lower redshift (2 < 1). 

(v) The SMF of galaxies alone is not sufficient for con¬ 
straining the galaxy formation models and the M* — /mol 
relation forms one of the important observational constraints 
on feedback processes implemented in galaxy-formation 
models. 

It is to be noted that we adopted the indirect estimate s of 
/mol using the optical survey data llPopping et al.|[2012ll for 
comparison with the galaxy formation models in this study. 
To form a definitive conclusion regarding the /mol evolution 
of galaxies (particularly less massive and massive galaxies), 
it is necessary to directly measure the cold gas mass of galax¬ 
ies over wide ranges of stellar mass and redshift. The mass 
function of the stellar component is of course a basic ob¬ 
served constraint, which must be satisfied for the models. 
The cold-gas component is also an important component in 
the baryonic cycle of galaxies. A combination of the obser¬ 
vational constraints on stellar and cold gas components will 
aid in sketching a complete scenario of galaxy formation and 
evolution. Observations with forthcoming facilities such as 
the ALMA, the NOEMA, and the SKA, will allow us to 
explore the statistical properties and evolution of cold-gas 
components in galaxies, such as Hi and H 2 mass functions, 
as well as the M* — /mol relation. 
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(a) z<0.1 (b) 0.1<z<0.5 (c) 1.0<z<1.5 (d) 2.0<z 



Figure 4. Total gas fraction /gas as a function of stellar mass and redsliift (a: z < 0.1, b: 0.1 < z < 0.5, c: 1.0 < z < 1.5 and d: 
2.0 < z). The black a nd grey dot-dash lines indicate the M* — /gas relation based on the K-S law for whole and star-forming samples in 
iPoDping et al.l ll2012h . The redshifts of the equations for /mol,Pl2,SF S'lid /inol,Pl2 3 -re 0. 025 for (a), 0.3 for (b), 1.25 for (c), and 2.2 for (d). 
The lines with filled triangles indicate the results of the ‘equilibrium model (EM)’ bv IPave et al.l ll2012h (red solid lines: model with all 
feedbacks, light-blue dot-dashed lines: model without photoionising feedback (PHOTO), green dotted lines: model without shock-heating 
during accretion (GRAY), purple short-dotted lines: model without quenching associated with black hole growth (QUENCH), and blue 
dashed lines: model without additional feedback associated with wind (WIND)). On each plot, we show the result of EM calculations at 
z = 0 in (a), z = 0.2 in (b), z = 1 in (c) and z = 2 in (d). The Mj, — /gas re l ations from the direc t observations (HI and CO) are plotted 
as a reference llLerov et al.ll2003 : ISaintonge et al.ll201l3 : iBoselli et al.ir2014l : iBothwell et al.ll2014ll . 


productive discussions, and Ross Burns, Daisuke lono and 
Editage (www.editage.jp) for English language editing. 

Funding for SDSS-III has been provided by the Al¬ 
fred P. Sloan Foundation, the Participating Institutions, 
the National Science Foundation, and the U.S. Depart¬ 
ment of Energy Office of Science. The SDSS-III web site 
is (http://www.sdss3.org/). SDSS-III is man- aged by the 
Astrophysical Research Consortium for the Participating 
Institutions of the SDSS-III Collaboration including the 
University of Arizona, the Brazilian Participation Croup, 
Brookhaven National Laboratory, Carnegie Mellon Uni¬ 
versity, University of Florida, the French Participation 
Group, the German Participation Group, Harvard Univer¬ 


sity, the Institute de Astrofisica de Canarias, the Michigan 
State/Notre Dame/JINA Participation Group, Johns Hop¬ 
kins University, Lawrence Berkeley National Laboratory, 
Max Planck Institute for Astrophysics, Max Planck Insti¬ 
tute for Extraterrestrial Physics, New Mexico State Univer¬ 
sity, New York University, Ohio State University, Pennsyl¬ 
vania State University, University of Portsmouth, Princeton 
University, the Spanish Participation Group, University of 
Tokyo, University of Utah, Vanderbilt University, University 
of Virginia, University of Washington, and Yale University. 
































Redshift evolution of M* — /moi relation int)<z<2 13 


REFERENCES 

Abazajian K. N. et al., 2009, ApJS, 182, 543 
Ahn C. P. et al., 2014, ApJS, 211, 17 
Arnett D., 1996, Supernovae and Nucleosynthesis: An In¬ 
vestigation of the History of Matter from the Big Bang to 
the Present 

Bauermeister A., Blitz L., Bolatto A., Bureau M., Teuben 
P., Wong T., Wright M., 2013, ApJ, 763, 64 
Baugh C. M., Lacey C. G., Frenk C. S., Granato G. L., Silva 
L., Bressan A., Benson A. J., Gole S., 2005, MNRAS, 356, 
1191 

Benson A. J., Bower R. G., Frenk G. S., Lacey G. G., Baugh 
G. M., Cole S., 2003, ApJ, 599, 38 
Berry M., Somerville R. S., Haas M. R., Gawiser E., Mailer 
A., Popping G., Trager S. C., 2014, MNRAS, 441, 939 
Bigiel F., Leroy A., Walter F., Brinks E., de Blok W. J. G., 
Madore B., Thornley M. D., 2008, AJ, 136, 2846 
Bird S., Vogelsberger M., Haehnelt M., Sijacki D., Genel 
S., Torrey P., Springel V., Hernquist L., 2014, MNRAS, 
445, 2313 

Blanton M. R., Roweis S., 2007, AJ, 133, 734 
Blitz L., Rosolowsky E., 2004, ApJL, 612, L29 
Blitz L., Rosolowsky E., 2006, ApJ, 650, 933 
Bolatto A. D., Wolfire M., Leroy A. K., 2013, ARA&A, 51, 
207 

Booth C. M., Schaye J., 2009, MNRAS, 398, 53 
Boselli A., Cortese L., Boquien M., 2014, A&A, 564, A65 
Bothwell M. S. et al., 2014, MNRAS, 445, 2599 
Bouche N. et ah, 2010, ApJ, 718, 1001 
Bower R. G., Benson A. J., Malbon R., Helly J. C., Frenk 
C. S., Baugh C. M., Cole S., Lacey C. G., 2006, MNRAS, 
370, 645 

Brinchmann J., Chariot S., White S. D. M., Tremonti C., 
Kauffmann G., Heckman T., Brinkmann J., 2004, MN¬ 
RAS, 351, 1151 

Bruzual G., Chariot S., 2003, MNRAS, 344, 1000 
Carilli C. L., Walter F., 2013, ARA&A, 51, 105 
Catinella B., Haynes M. P., Giovanelli R., Gardner J. P., 
Connolly A. J., 2008, ApJL, 685, L13 
Chabrier G., 2003, PASP, 115, 763 
Chomiuk L., Povich M. S., 2011, AJ, 142, 197 
Christensen C., Quinn T., Governato F., Stilp A., Shen S., 
Wadsley J., 2012, MNRAS, 425, 3058 
Ciesla L. et ah, 2014, A&A, 565, A128 
Cole S., Lacey C. G., Baugh C. M., Frenk C. S., 2000, 
MNRAS, 319, 168 

Cole S. et al., 2001, MNRAS, 326, 255 
Conroy C., Wechsler R. H., 2009, ApJ, 696, 620 
Cortese L. et ah, 2012, A&A, 544, AlOl 
Cortese L. et ah, 2014, MNRAS, 440, 942 
Croton D. J. et al., 2006, MNRAS, 365, 11 
Daddi E. et al., 2010, ApJ, 713, 686 

Daddi E., Cimatti A., Renzini A., Fontana A., Mignoli M., 
Pozzetti L., Tozzi P., Zamorani G., 2004, ApJ, 617, 746 
Dalla Vecchia C., Schaye J., 2008, MNRAS, 387, 1431 
Dalla Vecchia C., Schaye J., 2012, MNRAS, 426, 140 
Dannerbauer H., Daddi E., Riechers D. A., Walter F., Car¬ 
illi C. L., Dickinson M., Elbaz D., Morrison G. E., 2009, 
ApJL, 698, L178 

Dave R., Finlator K., Oppenheimer B. D., 2011a, MNRAS, 
416, 1354 


Dave R., Finlator K., Oppenheimer B. D., 2012, MNRAS, 
421, 98 

Dave R., Oppenheimer B. D., Finlator K., 2011b, MNRAS, 
415, 11 

De Lucia G., Blaizot J., 2007, MNRAS, 375, 2 
Dekel A., Sari R., Ceverino D., 2009, ApJ, 703, 785 
Dekel A., Silk J., 1986, ApJ, 303, 39 
Duffy A. R., Kay S. T., Battye R. A., Booth C. M., Dalla 
Vecchia C., Schaye J., 2012, MNRAS, 420, 2799 
Efstathiou G., 1992, MNRAS, 256, 43P 
Fakhouri O., Ma C.-P., Boylan-Kolchin M., 2010, MNRAS, 
406, 2267 

Faucher-Giguere C.-A., Keres D., Ma C.-P., 2011, MNRAS, 
417, 2982 

Fernandez X. et al., 2013, ApJL, 770, L29 
Forster Schreiber N. M. et al., 2009, ApJ, 706, 1364 
Freudling W. et al., 2011, ApJ, 727, 40 
Fu J., Guo Q., Kauffmann G., Krumholz M. R., 2010, MN¬ 
RAS, 409, 515 

Fu J., Kauffmann G., Li C., Guo Q., 2012, MNRAS, 424, 
2701 

Fumagalli M., Krumholz M. R., Hunt L. K., 2010, ApJ, 
722, 919 

Gabor J. M., Dave R., Oppenheimer B. D., Finlator K., 
2011, MNRAS, 417, 2676 

Geach J. E., Smail L, Moran S. M., MacArthur L. A., Lagos 

C. d. P., Edge A. C., 2011, ApJL, 730, L19 

Geach J. E., Smail L, Moran S. M., Treu T., Ellis R. S., 
2009, ApJ, 691, 783 

Genel S. et al., 2014, MNRAS, 445, 175 
Genzel R. et ah, 2010, MNRAS, 407, 2091 
Gnedin N. Y., 2000, ApJ, 542, 535 
Gnedin N. Y., Kravtsov A. V., 2011, ApJ, 728, 88 
Gnedin N. Y., Tassis K., Kravtsov A. V., 2009, ApJ, 697, 
55 

Haardt F., Madau P., 2001, in Clusters of Galaxies and 
the High Redshift Universe Observed in X-rays, Neumann 

D. M., Tran J. T. V., eds., p. 64 

Heifer T. T., Thornley M. D., Regan M. W., Wong T., 
Sheth K., Vogel S. N., Blitz L., Bock D. C.-J., 2003, ApJS, 
145, 259 

Henriques B. M. B., White S. D. M., Thomas P. A., Angulo 
R., Guo Q., Lemson G., Springel V., Overzier R., 2015, 
MNRAS, 451, 2663 

Henriques B. M. B., White S. D. M., Thomas P. A., Angulo 
R. E., Guo Q., Lemson G., Springel V., 2013, MNRAS, 
431, 3373 

Ikeuchi S., 1986, Ap&SS, 118, 509 
Ilbert O. et ah, 2010, ApJ, 709, 644 
Kauffmann G., 1996, MNRAS, 281, 475 
Kauffmann G., Colberg J. M., Diaferio A., White S. D. M., 
1999, MNRAS, 303, 188 
Kauffmann G. et ah, 2003, MNRAS, 346, 1055 
Kauffmann G. et ah, 2012, MNRAS, 422, 997 
Kennicutt R. C., Evans N. J., 2012, ARA&A, 50, 531 
Kennicutt, Jr. R. C., 1989, ApJ, 344, 685 
Kennicutt, Jr. R. C., 1998a, ARA&A, 36, 189 
Kennicutt, Jr. R. C., 1998b, ApJ, 498, 541 
Keres D., Katz N., Weinberg D. H., Dave R., 2005, MN¬ 
RAS, 363, 2 

Khochfar S., Ostriker J. P., 2008, ApJ, 680, 54 
Kochanek C. S. et ah, 2001, ApJ, 560, 566 


14 K. Morokuma-Matsui and J. Baba 


Kravtsov A. V., Gnedin O. Y., Klypin A. A., 2004, ApJ, 
609, 482 

Kroupa R, 2001, MNRAS, 322, 231 
Krumholz M. R., 2013, MNRAS, 436, 2747 
Krumholz M. R., McKee C. F., Tumlinson J., 2008, ApJ, 
689, 865 

Krumholz M. R., McKee C. F., Tnmlinson J., 2009a, ApJ, 
693, 216 

Krumholz M. R., McKee C. F., Tumlinson J., 2009b, ApJ, 
699, 850 

Lacey C., Cole S., 1993, MNRAS, 262, 627 
Lagos C. D. P., Baugh C. M., Lacey C. G., Benson A. J., 
Kim H.-S., Power C., 2011a, MNRAS, 418, 1649 
Lagos C. d. P. et al., 2015, ArXiv e-prints 
Lagos C. D. P., Lacey C. G., Baugh C. M., Bower R. G., 
Benson A. J., 2011b, MNRAS, 416, 1566 
Larson R. B., 1974, MNRAS, 169, 229 
Leitner S. N., 2012, ApJ, 745, 149 

Leroy A. K., Walter F., Brinks E., Bigiel F., de Blok 
W. J. G., Madore B., Thornley M. D., 2008, AJ, 136, 
2782 

Lu Y., Mo H. J., Wechsler R. H., 2015, MNRAS, 446, 1907 
Madau P., Dickinson M., 2014, ARA&A, 52, 415 
Magdis G. E. et ah, 2012a, ApJ, 760, 6 
Magdis G. E. et ah, 2012b, ApJL, 758, L9 
Makiya R., Totani T., Kobayashi M. A. R., Nagashima M., 
Takeuchi T. T., 2014, MNRAS, 441, 63 
Mancini C. et ah, 2011, ApJ, 743, 86 
Martig M., Bournaud F., Teyssier R., Dekel A., 2009, ApJ, 
707, 250 

Martig M. et ah, 2013, MNRAS, 432, 1914 
Mitchell P. D., Lacey C. G., Cole S., Baugh C. M., 2014, 
MNRAS, 444, 2637 

Mo H. J., Yang X., van den Bosch F. C., Katz N., 2005, 
MNRAS, 363, 1155 

Moran S. M., Ellis R. S., Treu T., Smith G. P., Rich R. M., 
Smail L, 2007, ApJ, 671, 1503 
Morokuma-Matsui K., Baba J., Sorai K., Kuno N., 2015, 
PASJ, 67, 36 

Nagashima M., Yahagi H., Enoki M., Yoshii Y., Gouda N., 
2005, ApJ, 634, 26 

Narayanan D., Bothwell M., Dave R., 2012a, MNRAS, 426, 
1178 

Narayanan D., Krumholz M. R., Ostriker E. C., Hernquist 
L., 2012b, MNRAS, 421, 3127 
Navarro J. F., Steinmetz M., 1997, ApJ, 478, 13 
Noeske K. G. et ah, 2007, ApJL, 660, L43 
Obreschkow D., Croton D., De Lucia G., Khochfar S., 
Rawlings S., 2009, ApJ, 698, 1467 

Okamoto T., Gao L., Theuns T., 2008, MNRAS, 390, 920 
Okamoto T., Shimizu I., Yoshida N., 2014, PASJ, 66, 70 
Oppenheimer B. D., Dave R., Keres D., Fardal M., Katz 
N., Kollmeier J. A., Weinberg D. H., 2010, MNRAS, 406, 
2325 

Popping G., Behroozi P. S., Peeples M. S., 2015, MNRAS, 
449, 477 

Popping G., Caputi K. L, Somerville R. S., Trager S. C., 
2012, MNRAS, 425, 2386 

Popping G., Somerville R. S., Trager S. C., 2014, MNRAS, 
442, 2398 

Rahmati A., Pawlik A. H., Raicevic M., Schaye J., 2013a, 
MNRAS, 430, 2427 


Rahmati A., Schaye J., 2014, MNRAS, 438, 529 
Rahmati A., Schaye J., Bower R. G., Crain R. A., Furlong 
M., Schaller M., Theuns T., 2015, ArXiv e-prints 
Rahmati A., Schaye J., Pawlik A. H., Raicevic M., 2013b, 
MNRAS, 431, 2261 

Rees M. J., Ostriker J. P., 1977, MNRAS, 179, 541 
Rhee J., Zwaan M. A., Briggs F. H., Chengalur J. N., Lah 
P., Oosterloo T., van der Hulst T., 2013, MNRAS, 435, 
2693 

Saintonge A. et al., 2011a, MNRAS, 415, 32 
Saintonge A. et al., 2011b, MNRAS, 415, 61 
Saintonge A. et al., 2013, ApJ, 778, 2 
Santini P. et ah, 2014, A&A, 562, A30 
Schaye J. et al., 2015, MNRAS, 446, 521 
Schaye J., Dalla Vecchia C., 2008, MNRAS, 383, 1210 
Schaye J. et al, 2010, MNRAS, 402, 1536 
Scoville N. et al., 2007, ApJS, 172, 1 
Scoville N. et al., 2014, ApJ, 783, 84 
Scoville N. et al., 2015, ArXiv e-prints 
Sijacki D., Springel V., Di Matteo T., Hernquist L., 2007, 
MNRAS, 380, 877 

Silk J., Mamon G. A., 2012, Research in Astronomy and 
Astrophysics, 12, 917 

Somerville R. S. et al., 2008a, ApJ, 672, 776 
Somerville R. S., Dave R., 2014, ArXiv e-prints 
Somerville R. S., Gilmore R. C., Primack J. R., Dominguez 
A., 2012, MNRAS, 423, 1992 

Somerville R. S., Hopkins P. F., Cox T. J., Robertson B. E., 
Hernquist L., 2008b, MNRAS, 391, 481 
Somerville R. S., Kolatt T. S., 1999, MNRAS, 305, 1 
Somerville R. S., Popping G., Trager S. C., 2015, ArXiv 
e-prints 

Somerville R. S., Primack J. R., 1999, MNRAS, 310, 1087 
Springel V., 2010, MNRAS, 401, 791 
Springel V., Di Matteo T., Hernquist L., 2005a, MNRAS, 
361, 776 

Springel V., Hernquist L., 2003, MNRAS, 339, 289 
Springel V. et ah, 2008, MNRAS, 391, 1685 
Springel V. et ah, 2005b, Nature, 435, 629 
Sutherland R. S., Dopita M. A., 1993, ApJS, 88, 253 
Tacconi L. J. et ah, 2010, Nature, 463, 781 
Tacconi L. J. et ah, 2013, ApJ, 768, 74 
Thompson R., Nagamine K., Jaacks J., Choi J.-H., 2014, 
ApJ, 780, 145 

Tomassetti M., Porciani C., Romano-Diaz E., Ludlow 
A. D., 2015, MNRAS, 446, 3330 
Tremonti C. A. et al., 2004, ApJ, 613, 898 
Vogelsberger M., Genel S., Sijacki D., Torrey P., Springel 
V., Hernquist L., 2013, MNRAS, 436, 3031 
Vogelsberger M. et al., 2014, MNRAS, 444, 1518 
White C. E., Somerville R. S., Ferguson H. C., 2015, ApJ, 
799, 201 

Wiersma R. P. C., Schaye J., Smith B. D., 2009, MNRAS, 
393, 99 

Wuyts S. et ah, 2011, ApJ, 738, 106 

Zibetti S., Chariot S., Rix H.-W., 2009, MNRAS, 400, 1181 

This paper has been typeset from a Tj^X/ DTj^X file prepared 
by the author. 


